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Missing energy signals such as monojets are a possible signature of Dark Matter (DM) at collid¬ 
ers. However, neutrino interactions beyond the Standard Model may also produce missing energy 
signals. In order to conclude that new “missing particles” are observed the hypothesis of BSM neu¬ 
trino interactions must be rejected. In this paper, we first derive new limits on these Non-Standard 
neutrino Interactions (NSIs) from LHC monojet data. For heavy NSI mediators, these limits are 
much stronger than those coming from traditional low-energy v scattering or v oscillation experi¬ 
ments for some flavor structures. Monojet data alone can be used to infer the mass of the “missing 
particle” from the shape of the missing energy distribution. In particular, 13 TeV LHC data will 
have sensitivity to DM masses greater than ~ 1 TeV. In addition to the monojet channel, NSI can 
be probed in multi-lepton searches which we find to yield stronger limits at heavy mediator masses. 
The sensitivity offered by these multi-lepton channels provide a method to reject or confirm the DM 
hypothesis in missing energy searches. 


I. INTRODUCTION 

Missing energy signals are the tell-tale clue of the pro¬ 
duction of stable neutral objects. Indeed the imbalance 
of momentum and energy is in fact precisely the way 
in which the neutrino was first discovered. Supposing 
that the LHC finds anomalous “missing energy” events 
above SM backgrounds, the determination of its origin 
will be of paramount importance. As known sources of 
missing energy, a plausible origin of new missing energy 
data will be neutrinos. However, new neutral particles 
beyond Standard Model (BSM) such as dark matter can 
also produce missing energy signals at colliders. 

In this paper we explore how LHC data can be used 
to distinguish these two potential sources of missing en¬ 
ergy. We illustrate that both the DM mass and the SU (2) 
charge of neutrinos can be used to discriminate between 
singlet DM and SM neutrinos. ^ For simplicity we will 
focus on the so-called “monojet” signature in which a 
single hard jet recoils against “nothing” mso]. Previous 
work in this direction using Tevatron and early LHC data 
was carried out in m- 

We begin by reviewing current experimental limits on 
neutrino-proton interactions in the context of effective 
field theory (EFT). Up to dimension 6, we can have 

• Neutrino magnetic dipole moments: 

if D (1) 


^ If DM itself transforms non-trivially under SU(2) the situation 
is more complex. We leave for future work a systematic study 
in this direction but note that some of the implications of SU{2) 
charged DM in a variety of representations has been studied in 
e.g. [T]. 


where the spin matrix is = i /2, Uv is 

the magnetic moment (measured in units of the 
Bohr magneton fiB = e/ (2me), where e, me are 
the charge and mass of the electron). 

• Non-standard neutrino interactions (NSIs), 

^NSi = -2V2 GfsI^^ {VaApyp) ijYPf) ■ ( 2 ) 

f P 

where the matrix specifies the strength of the 
v-f interaction, in units of Fermi’s constant, Gp = 
\/y/2v%w ~ 1.2 X 10-^ GeV■^ with vew = 246 
GeV. The labels a, j3 are flavor indices running over 
e, pL, T, and P is a projection operator. We take 
/ to be any SM fermion (though only the vector 
components oi f = e,u,d are relevant for neutrino 
oscillations). 

Let us first consider whether neutrino magnetic mo¬ 
ments below currents limits can produce sizeable missing 
energy at the LHC? For Majorana neutrinos the 3x3 
matrix does not have diagonal entries and is anti¬ 
symmetric, but is completely general if they are instead 
Dirac. In the SM the magnetic moment is proportional 
to the neutrino mass, and therefore extremely small, 
^SM ^ 10“20 Pqj- piraQ neutrinos in BSM scenarios, 
naturalness considerations on the coeffecients of effective 
operators imply, [31], far below present 

experimental sensitivity. Finally for Majorana neutrinos 
reactor data as measured by the GEMMA spectrometer 
constrains, fXi, < 3.2 x 10“^^ /xb, |32| while the 7 TeV 
LHC sensitivity is around ~ 3 x 10“^ /xb [33], far above 
what is allowed by reactor and solar data [31]. We con¬ 
clude that neutrino magnetic moments will not produce 
sizeable missing energy at the LHC. 

Proceeding now to operators of mass dimension 6, we 
turn our attention to the NSI operators between quarks- 
neutrinos. Non-standard neutrino interactions (NSIs) 
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FIG. 1: Typical Feynman diagrams for pp — >■ vv{XX) + j {left panel) and pp ^ vu ^ + v {right panel). Though 

singlet DM and neutrinos are largely degenerate in the former process, only SM neutrinos give rise to the latter process. 


were first introduced in 1977 [3S] and continue to be of 
wide phenomenological interest [11] [36H48] (see |36l 09] 
for reviews). 

They are constrained by solar [371 [5Mn| , atmo¬ 
spheric 23 IMl mi [STHSO] , long-baseline |3TH33l 05] 03 
IbUl [61], collider [TT] 04] 071 [62], cosmological [63], and 
neutrino scattering data 23 00100) • 

The Lorentz structure of Eq. [^ can be understood as 
follows. First, assume that NSI can be parameterized 
as Onsi = where are neutrino and 

SM fermion bilnears. Under the assumption that lepton 
number remains a good symmetry and only left-handed 
neutrinos enter into Gnsi-. all such operators can be de¬ 
composed into {V — A) ® {V — A), {V — A) ® {V + A) 
components. 

One may worry that sizeable NSI would also induce 
large charged lepton interactions 23 IMl IS5]- Indeed, 
to evade the very strong limits from the charged lepton 
equivalent of Eq. ([^ we consider dimension-8 operators 
of the form [62] 


r^dim—8 
°^NSI 


4e' 


fP 

Ot^ 


'^EW 


(HLaJtiHLfi) {q-ff,q) 


(3) 


where H is the SM Higgs doublet. In unitary gauge and 
upon electroweak (EW) symmetry breaking we can make 
the replacement H ^ {h + vew) /v^- Thus at low ener¬ 
gies, one indeed generates Eq. ([^ without charged lepton 
interactions of the same strength. 

The remainder of this paper is organized as follows. 
First we introduce our simplified model and calculational 
framework in Sec. im In Section IIIII we derive new con¬ 
straints on NSI based on the latest monojet data from 
the LHC. Then we turn to projections of monojet sensi¬ 
tivity at 13 TeV and the ability to infer the mass of the 
’’missing particles” from the shape of the distribu¬ 
tion. We find that for contact interactions, DM masses 
> 700 GeV can be discriminated from NSI with about 
100 fb“^ of integrated luminosity. In Sect. IV we then use 


two distinct multi-lepton channels to probe NSI. These 
searches have neutrino flavor dependent sensitivity and 
have better sensitivity than mono jets for heavy media¬ 
tors of NSI. In Sec. [V]we discuss the complementarity of 


these channels along with low-energy probes of NSI for 
DM-neutrino discrimination and conclude in Sec. I VII 


II. MODEL AND CALCULATIONAL 
FRAMEWORK 

In order to derive LHC limits on NSI/DM couplings e 
we have implemented two models in the Universal Feyn- 
Rules Format (UFO) [66] by adding to the SM a spin-1 
mediator, , which interacts with neutrinos, quarks and 
DM X through the phenomenological Lagrangians: 


■^NSI = 9u {vPLjf,v) + {qinidg + R^, 

= 9x RP + {qji_c{9q + 9^^^)^) R^^ 

+ mxXX (4) 

where v and q are summed over all neutrino and quark 
flavors respectively, and mx is the DM mass. The La- 
grangian Afjqsi correctly reproduces the contact interac¬ 
tion, Eq. ([^ when the vector mass, mu is large compared 
to the center of mass energy. Note that the DM literature 
tends to report limits on the scale of the dimension-six 
operator. A, defined as {X^f^X){q"fPq)/. The conver¬ 
sion from A to the NSI e parameter in this context is, 
e = (2GFA2)-b 

The main aim of this paper is to illustrate how 
can be discriminated from and gauge the relevant 

parametric dependencies present in s-channel comple¬ 
tions of NSI (t-channel completions are very strongly con¬ 
strained [1307] and not considered further). For details 
on a more complete Z' model we refer the reader to the 
Appendix and to [HTIIM] for additional models. Further¬ 
more, it is important to highlight that a complete model 
typically produces signatures in addition to the monojet 
and multilepton channels we consider here, making our 
approach conservative. 

Simplified models of the type in Eq. [^ have been stud¬ 
ied extensively in the DM literature [T] |3113 HU 113 1771 
IMl I 6 M 73 . While dijet searches provide additional con¬ 
straints on the models considered here (see e.g. [I3l IM] 1 
both .i^NSi and -SfoM contribute equally to this channel 
and thus it is not a useful discriminatory tool. 
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FIG. 2: Here we display the CMS monojet limits [88] on NSI 
at 95%CL for three different choices of the mediator width at 
= 8TeV and with integrated luminosity £ = 19.5 fb~^. 
The black solid line denotes the expected limit at 95%CL with 
y^s = 13TeV and L — 20fb“^. 


Our calculational framework is as follows: To keep the 
analysis simple, we consider only vector couplings, i.e. 
dt ~ 9x ~ 0. We import the UFO model into the 
MadGraph5_aMC@NLD framework m, where helicity am¬ 
plitudes are generated by the ALOHA (81] code. The hard 
scattering simulation is then processed through parton 
showering and hadronization using Pythia 6 |82] and 
Pythia 8 |13|. Finally we performed a fast detector sim¬ 
ulation for the monojet analysis, using both the PGS [84] 
and DELPHES |5S| programs to check our results. 

Eor the monojet computation we use the CTEQ 6L1 
[5S] set of parton distribution functions, as this is 
used by the experimental collaboration, and NNPDF 
2.3 |57j for the other processes. We chose the default 
dynamical factorization and renormalization scales of 
MadGraph_aMC(§NLO. 


III. MONOJET SEARCHES 

Any long-lived or stable neutral states, such as neu¬ 
trinos and DM, with couplings to protons can lead to 
monojet events at the LHC. These monojet processes, 
depicted in Fig.(le/t), are characterized by large miss¬ 
ing transverse energy and a very hard jet. In | 88 | . the 
CMS experiment searched for monojets with y^s = 8 TeV 
in the center of mass energy and C — 19.5 fb~^ of inte¬ 
grated luminosity, reporting an upper limit at 90%CL of 
e = 0.053 for a vector operator and an invisible particle 
mass mx = 1 CeV. 

To estimate the NSI signal we compute the cross sec¬ 


tions for the hard scattering process 

pp —i' V —^Vu + 1, 2j, (5) 

with one and two jets (quarks or gluons), j using the 
framework described in sec. ©• In particular we use 
the MLM prescription [521 for matching matrix elements 
with soft jets from the parton shower. Following the CMS 
analysis [HH| we require the leading jet to have prij) > 
110 CeV and to be in the central region of the detector 
|? 7 (j)| < 2.6. Events with more than 3 jets with px > 
30 CeV and |r?| < 4.5 are discarded, while a second jet is 
allowed as long as the difference in azimuthal angle to the 
leading jet is less than 2.5, A0(ji, < 2.5. We further 
require the missing transverse energy 0x > 450 CeV, 
found to give the best discriminant. 

With this analysis set-up we found excellent agreement 
in the shape of the missing transverse energy distribution 
for Z{vv)+ jets and W{(.v)+ jets SM background. We 
also found agreement within scale and PDE uncertainties 
for the number of events. We nevertheless use the fact 
that the collaboration provides a more precise prediction 
from data driven techniques and we rescale our predic¬ 
tions by a correction factor of 1.19 to agree with their 
prediction. 

The CMS collaboration report 157 events as the upper 
95% confidence level (CL) limit on the number of events 
from new physics. Note that a downward fluctuation in 
the observed number of events gives a constraint about 
30% stronger than expected. We compute the resulting 
NSI limits that are shown in Eig. as a function of the 
mediator mass and width (Fy = ^). 

Note that flavor diagonal NSIs interfere with the dom¬ 
inant SM background process, pp ^ Z + j ^ Vu+ jets. 
The strength of the effect depends on the Lorentz struc¬ 
ture of the coupling, and the mass of the mediator. The 
effect is small in the contact interaction limit, < 5%, but 
can be as large as 20% when the mass of the mediator is 
close to the Z mass. Although interference is a feature 
specific to the NSI case it only affects the total num¬ 
ber of events and does not aid in distinguishing between 
dark matter and NSI. We shall therefore omit it in the 
following 

A. Projection to y^ = 13 TeV LHC and jet px 
shape analysis 

The next LHC run at i/s = 13 TeV will either further 
limit or discover NSI and/or DM in monojet searches. 
In Fig. we show our projected LHC monojet 95%CL 
limit as the solid black line for mx = 0 CeV at y/s = 
13 TeV with the luminosity £ = 20 fb“^ as expected for 
the first year of collisions. We use the same set-up used 
at ^/s = 8 TeV and the same normalization rescaling. 
We assume a systematic error of 5% and performed a 
analysis with which the expected bounds at 8 TeV 
quoted by CMS were reproduced within error. At this 
luminosity the systematic error dominates and increasing 
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missing Ej [GeV] 

FIG. 3: Distribution of events in missing transverse energy, 
^T, for ^ = 13TeV and C = 100 fb“^ for DM masses mx = 
0 GeV, 500 GeV, 700 GeV and 1 TeV. Here each distribution 
is generated assuming contact interactions which produce an 
identical total number of events, with an interaction strength 
just below present bounds (see e.g. Fig. [^. 


the luminosity further does not appreciably change the 
experimental sensitivity. 

The first observable we use to distinguish NSI from 
DM is the monojet 0t distribution. Sufficiently heavy 
DM masses are kinematically relevant at LHC energies 
and affect the shape of the distribution. 

The allowed value of e just below the 95%CL present 
limit is ~ 0.04 in the contact interaction limit, for 
massless missing energy particles. This situation, i.e. 
e = 0.04, would produce a 2.9 (t excess according to our 
projections for the first year Run II of the LHC. The 
same excess of events can be prodnced for lower e but 



FIG. 4: projection analysis. Here we generate events at 

yG = 13 TeV in a model with contact interactions just below 
present bounds with a 0 GeV DM mass. The 0 GeV] 

is computed by fitting mx = 500, 700,1000 GeV DM masses 
to the input data from a massless invisible particle. The 3 (t, 
4cr and 5a confidence levels are plotted for reference. 


lighter mediator mass or larger e and heavier particles in 
the final state. In Fig. we show the distribution 
for L = lOOfb”^ for DM masses, mx = 500,700,1000 
GeV. The total cross sections are all normalized to the 
e = 0.04 massless case (shown in the projection Fig.[^ so 
that all signals produce the same total number of events 
with JS't > 450 GeV. 

The shape of the distributions clearly allow to dis¬ 
tinguish between the heaviest and lightest DM masse. 
We quantify this in a simple analysis. In addition to 
the statistical error, we assume a systematic error per bin 
of 5% for the background and 20% for the dark matter 
contribution, as reported in [88]. The distribution is 
then given by 


X^[tox; 0 GeV] 


E 


S,(mx) - S,(0GeV) 
o-i 


( 6 ) 


where Si(rnx) is the number of events in the i**' bin. 
The distribution is shown in Fig. [^for the three masses 


mx = 500,700,1000. 


IV. MULTILEPTON SEARCHES 

In addition to the monojet signal, NSI can produce 
signals in other channels due to the SU{2) charge of neu¬ 
trinos. For example, as shown in Fig. [^one of the pro¬ 
duced neutrinos can radiate a W boson that decays to 
either jets or £ + v, 

pp ^ Vv ^ V + . (7) 

Mutli-lepton searches of this type have been used previ¬ 
ously to constrain NSI using LHC data mmi]. 

In order to exclude the NSI hypothesis and claim the 
discovery of a new source of missing energy, we must ex¬ 
clude all possible neutrino flavor structures of NSI. For 
this it is necessary to consider the lepton in the final state 
to be a tau, a muon or an electron. Since the mixed flavor 
interaction, e.g. Srfj,, £e/i, will regardless produce one of 
these leptons, this condition is also sufficient to constrain 
mixed terms. For the muon and electron in the final state 
we have relied on the ^/s = 8 TeV and £ = 20.3 fb~^ AT¬ 
LAS search for resonant diboson production where one 
boson decays leptonically and the other hadronically m- 
For the tau lepton final state we have used the ATLAS 
search for supersymmetry with large missing transverse 
energy, jets and at least one tau lepton, at ^/s = 8 TeV 
and £ = 20.3 fb“^ of data (SD]. We will briefly describe 
each analysis and results in the following. 

The searches we used are not optimized for the NSI 
signal topologies and we expect that dedicated analy¬ 
ses can improve our results. Moreover NSI can lead to 
signals not considered here but they are expected to be 
sub-dominant. For example pp vv^Z -x jj/t^£~), 
where the neutrino radiates a Z boson will suffer from 
large background from Drell-Yan production. Similarly, 
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1 + r + MET 


1 + W(-+jets) + MET, HRR 




FIG. 5: Each panel displaying individual LHC search limits on NSI for three different choices of the mediator width.The left 
panel displays the r + ^ + MET search from ATLAS [90], while the right panel shows the sensitivity from the jj + £ + MET 
search from ATLAS m (see text for details). 


in pp — i/£{W —>■ £v) with highly energetic £~^£~ sys¬ 
tem, the W cannot be reconstructed, suffering from many 
more backgrounds. Nonetheless, all these channels may 
contribute to put bounds on NSI and require a dedicated 
analysis. 

We begin by considering hadronic decays of the Ws, 
i.e 


A. pp ^ V + — >■ jj, £ = e, p 

In this analysis the TF-boson is required to be highly 
boosted to reduce hadronic backgrounds. Consequently 
the two jets from the W are likely to appear as a single jet 
making jet substructure techniques relevant. The parton- 
level computation was passed through parton-showering 
and hadronization using Pythia 8. 

We employ the event selection of the experimental 
analysis in |91| : Leptons are required to have transverse 
momentum pp > 25 GeV and |ry| < 2.5. Moreover, they 
are required to satisfy the following isolation criteria: 
the scalar sum of pp of tracks with > 1 GeV within 
AR = Arf + AefP' = 0.2 of the lepton track is required 
to be less than 15% of the lepton pp. The missing trans¬ 
verse energy, defined as the negative of the vectorial sum 
of the transverse momenta of all electrons, muons and 
jets within \r]\ < 4.9, is required to be Xt > 30 GeV. 

We cluster the jets with two different jet definitions 
provided by Fast jet 3.1.2 [52] • For the signal region 
where the W has large pp and the jets cluster into a 
single ’fat’ jet, we use the Gambridge algorithm. Other¬ 
wise, we use the anti-fc^ algorithm with R = 0.4. The 
ATLAS analysis [HI] defines three signal regions as: The 


merged region (MR), the high-pT resolved region (HRR) 
and the low-pr resolved region (LRR) respectively. In 
the MR the largest pp jet (J) is taken to represent a de¬ 
cayed IV-boson, if it fulfills pp{ J) > 400 GeV, |r 7 (J)| < 2 
and 65 GeV < m( J) < 105 GeV with the azimuthal angle 

difference between J and ^p satisfying A())(J,Xt) < 1- 
Additionally, the pp of the lepton and ^p system is re¬ 
quired to be pp{£Xt) > 400 GeV. 

If the event does not pass these cuts, we proceed to 
the resolved region, where the two leading 0.4 ani\-kp 
jets, ji ^2 reconstruct a decayed W boson if: |p(j)| < 2.8, 
65 GeV < m{jj) < 105 GeV and A<j){ji,^p) < 1. 
The HRR (LRR) is defined by pp{jj) > 300(100) GeV, 
pp{j) > 80(30) GeV and pp{£Xt) > 300(100) GeV. 

After normalizing with an approximate NLO K-factor 
of K=1.7 [93], we get reasonable agreement in all three 
regions for the number of events expected from the SM 
diboson background. Therefore we assume that this sim¬ 
ple analysis is accurate enough for our needs. 

We used the model described by Eq. @ to estimate 
the visible cross sections, as, and associated number of 
events, S = as £ ot the NSI signal, where the luminosity 
is £ = 20.3 fb“^. We rescale our prediction by a K- 
factor K = 1.2 to account for QGD corrections extracted 
from on-shell Z' production |94] . We moreover assumed 
a conservative flat theoretical error of 30% to account 
for PDF and scale uncertainty. The SM prediction for 
the total number of events and uncertainty, B A as, is 
161500 ± 2300, 870 ± 40 and 295 ± 22 for LRR, HRR 
and MR respectively and the observed number of events, 
= 157837, 801 and 295 respectively. We summed 
the errors in quadrature, ^ -I- (0.35')^ and 

estimate the 95% GL upper limit on S using a analysis. 










6 


solving for S the equation 

\ CTtOT / 

The resulting limits in terms of e are shown in Fig. 

Next we will consider leptonically decaying W bosons 

B. pp ^ u + W^T^, —>■ £u and pp ^ u + 

—>■ TV, T —hadrons, £ = e, p 

The signal region defined in the ATLAS search HO] 
relevant for our final state is referred as the r+lepton 
“GMSB signal” region, which requires a reconstructed 
hadronically decayed tau lepton and a single isolated elec¬ 
tron or muon. Non standard neutrino interactions involv¬ 
ing a tau lepton will contribute to this process, but other 
NSI flavour structures without a tau leptons will equally 
contribute when the W boson decays to a tau lepton and 
tau neutrino. 

In our analysis we assume the tau is reconstructed with 
70% of efficiency in this region, as reported in the analy¬ 
sis. In addition, we reproduce the kinematical cuts given 
therein: prii) > 25 GeV, pt{t) > 20GeV, lepton trans¬ 
verse mass, mT{t) > 100 GeV, defined by 

Wt(^) = \J‘^PT{^)y^T (l — cos(A(/)(^,^t))) (9) 

and meff > 1700 GeV, where 

Weff =Pt(^) +Pt(t)-I-Xt- ( 10 ) 

The 95% CL limit on the visible cross section provided 
by the ATLAS collaboration is 0.20 fb for r -|- e channel 
and 0.26 fb for the t + p channel. Using these numbers 
we find the 95%CL exclusion limit shown in Fig. as 
the blue dot-dashed line. The limit shown is for NSI 
involving a tau lepton, e^T, however the difference with 
respect to other flavour structures is small. For e^e it is 
only a few percent, and for it is about 15% due to the 
weaker experimental upper limit on the muon channel. 

V. DISCUSSION 

After deriving new limits on NSI, we finally asses what 
future LHC data can unveil. If anomalous missing energy 
events appear in the next run of the LHC, they will be 
consistent with either DM or NSI just at the border of 
the current constraints. If the events are due to TeV scale 
DM, then Xr shape analysis will be enough to rule out 
NSIs as the origin. If that is not the case, then we can still 
use multi-lepton channels to help discriminate between 
NSIs and DM. To illustrate this point we consider three 
distinct benchmark scenarios: 

• Benchmark A, (my,e) = (100 GeV,0.15). 

The LHC is not a particularly good environment 



FIG. 6: In the mediator mass-coupling plane (mv,e), we com¬ 
pare existing searches for NSI from neutrino-nucleus scatter¬ 
ing to the LHC mono-jet limits derived in this paper. The 
upper curve in the gray band depicts the current monojet 
limits, while the lower curve shows the 13 TeV projection 
with 100 fb“^. The dot-dashed curves represent the current 
multi-lepton constraints on NSI based on 8 TeV LHC data. 
Additional low-energy constraints on the NSI parameter ea 0 
include NuTeV’s constraint on |95| and CHARM’s con¬ 
straint on £ee [96]. For reference the constraint on Ett is suf¬ 
ficiently weak that it does not appear on the plot (see |36|L 


for discriminating neutrinos from DM in the light 
mediator limit. Although NuTeV’s constraint [95] 
on /r-flavored diagonal NSI shown as the orange 
dashed line in Fig. allows us to conclude that 
this particular flavor structure is not responsible for 
anomalous missing energy events, the other flavor 
structures have much weaker constraints and can¬ 
not be excluded as potential explanations of LHG 
monojet signals. NSI with r or e flavored interac¬ 
tions can simultaneously escape low-energy probes 
and multi-lepton searches at the LHC. Future data 
from dedicated low-energy experiments searching 
for Vf, — N OT Vr — N may help resolving this. 

• Benchmark B, {mv,s) = (1500 GeV, 8 x 10“^). 
Here, the LHC’s ability to discriminate NSI from 
the DM hypothesis is much more favorable given 
the strength of e, p flavored NSI limits. Thus al¬ 
though monojet data will be at the discovery level, 
r-flavored NSI is degenerate with a DM interpre¬ 
tation. However, since the r -|- e -I- MET search 
utilized in the present paper, is not optimized for 
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NSI it is possible that a dedicated analysis could 
help resolve this. 

• Benchmark C, {mv,s) = (5 TeV, 0.03). 

This final benchmark is the most optimistic, as 
the discrimination between NSI and DM is robust. 
This is because there are two sensitive probes of 
NSI since both the r + e + MET and jj + £ + MET 
channels yield stronger constraints than monojets. 
Thus for example, monojet data originating from 
this benchmark would already be excluded from be¬ 
ing of NSI origin with present multi-lepton data. 

Summarizing, the multi-lepton probes are crucial for 
distinguishing between DM and NSIs in benchmark C 
and partially in the case of benchmark B. For benchmark 
A input from additional low-energy experiments will be 
needed. These can be either DM or neutrino probes. 
For example, DM direct detection data can be used to 
determine the mass of the DM and bound the mediator 
mass [57ITO] . 

Alternatively, in the case of neutrinos constraints on 
neutrino scattering will improve shortly. Using the meth¬ 
ods outlined in m. the COHERENT |1()()1 llOlj collab¬ 
oration’s multi-target measurement of coherent elastic 
neutrino-nucleus scattering can be used to substantially 
strengthen the limits on NSI [100] from NuTeV [SS] and 
CHARM [96]. 

Finally, thanks to the modification of neutrino oscilla¬ 
tion probabilities that (vector) NSI induces, long-baseline 
and solar neutrino data will also further limit NSI. Future 
probes of NSI include long-baseline experiments such as 
NOz^A and DUNE |45|, as well as atmospheric data from 
IceCube DeepCore [16] , and solar neutrino data from DM 
direct detection experiments |102j . 

These complementary experimental searches will be 
tremendously useful in obtaining better sensitivity to 
NSI. 


VI. CONCLUSION 

If anomalous events with missing energy are found in 
the next run of the LHC, determining the nature of the 
“missing particles” will be of utmost importance. Given 
that neutrinos are the only confirmed source of miss¬ 
ing energy to date, a neutrino interpretation would be 
quite natural. Moreover, such non-standard neutrino in¬ 
teractions are rather weakly constrained and could well 
produce sizeable j -I- MET rates at the LHC. Here we 
investigated two useful tools that may aid in this dis¬ 
crimination: shape analysis of monojet data, and 

multi-lepton data. 

We found that NSI can be discriminated from DM 
based on shape analysis if the DM mass is > 1 TeV. 

Next, the SU{2) charge of neutrinos implies that NSI 
contributes in channels involving charged leptons. This 
gives a simple discriminant between neutrino explana¬ 
tions of missing energy from singlet DM. To this end 


we studied jj + MET and t -1- e -I- MET events to 
derive new limits on NSI. We have found that NSI me¬ 
diators with masses > 800 GeV can be fairly robustly 
discriminated from DM interpretations. This is because 
the above multi-lepton channels offer greater sensitivity 
at large mediator masses than monojets. In particular, 
for mediator masses greater than 1.5 TeV both channels 
are separately strong enough to discriminate between NSI 
and DM. Light mediator NSI remains hard to probe with 
LHG data, which underscores the importance of upcom¬ 
ing low-energy probes of NSI. 
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Appendix: A UV Model of NSI 


We would like to arrive at a simple completion of 
Eq. (§, which suggests a spin-1 completion given its 
Lorentz structure. Moreover, the main feature of Eq. @ 
is that it implies stronger quark-neutrino interactions 
than quark-charged lepton interactions. A simple way to 
achieve is through the “baryonic portal” |103j (though 
see also [104H107j L In this class of models the quarks 
are charged directly under a new 11(1)' gauge symme¬ 
try. In addition there are new SM singlet fermions which 
also carry nonzero 11(1)' charge. We will refer to these 
as “baryonic neutrinos” for simplicity. When ordinary 
and baryonic neutrinos mass mix, the SM neutrinos effec¬ 
tively inherit a small piece of the new interaction. Thus 
one needs 


U gvVfj, {Q^f,Q -f FyNBL-lXfiNvi (11) 


where V is the gauge boson of the U(l)' symmetry, N is 
a singlet fermion, vi is the baryonic neutrino, and (j) is 
a baryonic Higgs whose VEV provides a mass for the V. 
Grucially, once (j) develops a VEV it allows mass mixing 
with N and hence N, vi, and Va all mass mix. Note that 
in the above we have used the standard notation for the 

Lepton doublets, = i ), where a = e, y, r, 


Next, notice that when mv,mN are both large com¬ 
pared to the momentum flowing through the V propaga¬ 
tor we can integrate out both new states to write 


£2s/2Gf 



sin^( 26 »b), 


( 12 ) 


where 9}, is the mixing angle. 











An important theoretical constraint on the model 
comes from anomaly cancellation |108) . The least con¬ 
strained possibility is when the new fermions are vector¬ 
like under the SM gauge group, but chiral under U{1)b- 
Some of these fermion carry electric charge, meaning that 
they are very strongly constrained. The most conser¬ 
vative constraint comes from chargino searches at LEP 
and imply that these fermions be heavier than ^ 90 
GeV 

Since the gauge boson mass is mv = gv{4>) I and the 
vector-like fermions have a mass controlled by the new 
VEV, iTif = X{cj)), the lower limit on the mass of the 
fermion translates into an upper limit on the size of the 
gauge coupling (also assuming A < tt) 


9v 


2Xmv 

nif 


< 6.3 X 10"^ 


100 GeV 

ruf 



my \ 

1 GeV/ 


(13) 


Thus the anomaly considerations in Eq. 
constrain the NSI coupling as: 


(131 directly 


e < 310 sin2(26lt,) (14) 


Thus anomaly considerations on their own, are not a sig¬ 
nificant constraint on NSI models of this sort. However 
the LHG constraints found in this paper, £ < 10“^ —10“^, 
can be translated to a constraint on the mixing between 
the baryonic neutrino and a SM neutrino, sin(20{,) < 
(2-6) X 10-3. 
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